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I. INTRODUCTION 



Exactly solvable potentials are extensively applied 
in the investigation of many physical systems in con- 
densed matter, nuclear physics, quantum optics and 
solid state physics. Typical examples within this large 
class of potentials in non-relativistic quantum me- 
chanics include the Morse potential the Poschl- 
Teller potential and the Coulomb potential. In par- 
ticular, the Morse potential has provided a reason- 
ably suitable approximation for studying anharmonic- 
ity and bond dissociation of diatomic molecules at 
low vibrational levels. Recently, coherent states ||], 
q-deformation |j and the phase space and trans- 
port properties of coupled Morse potentials ||] have 
also been investigated. The Poschl- Teller potential on 
the other hand has always been explored for its short- 
range properties. Its generalized coherent states 
nonlinear properties (7|^ and supersymmetric exten- 
sion 1^ have also been extensively studied. In re- 
cent years, under the guise of supersymmetric quan- 
tum mechanics, the idea of shaped invariant potentials 
have also been mooted [|4-18|. In particular, a large 
class of such potential is the Natanzon class . 

Algebraic methods |2^,^ exploiting the underly- 
ing Lie symmetry and its associated algebras have 
been widely used to study many of these exactly 
solvable potentials, for instance, among many oth- 
ers, Darboux transformation, Infeld-HuU factoriza- 
tion, Mielnick factorization, SUSY quantum mechan- 
ics, inverse scattering theory [0 and intertwining 
technique Moreover, by exploiting an underly- 



ing g-deformed quantum algebraic symmetry, Spiri- 
donov has shown how a finite difference differ- 
ential equation can generate in the limiting cases a 
set of exactly solvable potentials. Recently, algebraic 
method based on nonlinear algebras has also 

been used to analyze these solvable potentials. The 
main purpose of this paper is to develop a new uni- 
fied approach to obtain transition operators for a large 
class of exactly solvable potentials. 

In section ||, we briefly sketch the main ideas be- 
hind the diagonalization process. We then sketch the 
main idea behind our method of getting shift opera- 
tors from Hamiltonians with arbitrary solvable poten- 
tials in section [II. In section [V, we illustrate our 



method by showing how most of the exactly solvable 
potentials can be obtained through the method. Fi- 
nally in section we summarize our technique and 
discuss some subtle points associated with it. 



II. METHOD 

Suppose we have a physical system with Hamilto- 
nian H , generalized coordinate operator, Q, and gen- 
eralized momentum operator, P, satisfying the com- 
mutation relations 

[H,Q] = QQ^{H,h,c,)+PYl^{H,h,Ci), (la) 
\H,P] = QQ2{H,h,c,)+PIi2{H,h,Ci), (lb) 

where li, i € S for some arbitrary index set, S, are 
some invariants of motion of the system and Ci, i € S 
are some constants. Denoting 
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{[H, Q], [H, P]) := [H, {Q, P)], 



(2) 



we find that we can recast the previous commutation 
equations more succinctly as 



{[H, Q], [H, P])U^[H, (Q, P)U], 



(3) 



where ?7 is a 2 x 2 matrix so that the commutation 
relations in eq(0) can be rewritten as 

[H, {Q, P)] = (Q, P) ®| ) ^ (Q, P)M. (4) 

By diagonalizing the matrix, M, on the right hand 
side of eq(^) , we see that 

[H, {Q, P)U] = (Q, P)UU-^MU 



oscillator is the simplest example. In general, it is pos- 
sible to replace the coordinate operator x and momen- 
tum operator p by generalized "coordinate" operator 
Q and "momentum" operator P which are functions 
of the operators p, H and £ S. Physically, the 
energy spectrum will no longer be uniformly spaced 
since the eigenvalues, Ei and E2, associated with the 
raising and lowering operators are in functions of H 
and the amount through which these operators raise 
or lower depend very much on the energy level as well. 
This is the main source of nonlinearity in the algebra. 
Moreover, we note that the coefficient matrix, M, in 
eq(^ as well as the corresponding transformation ma- 
trix, U , are generally functions of the Hamiltonian H , 
the invariants of motion /,- and some constants c, . 



givmg 



[H, (Q, P)U] - (g, P)U 



El 
E2 



where Ei and E2 are the eigenvalues of M . By defin- 
ing (Q, P)U as 82), we see that the operators Si 
and 5*2 acts as lowering and raising operators. 

In the specific case of the one dimensional har- 
monic oscillator the Hamiltonian is given by 
H = \{x'^ where x is the coordinate operator 
and p — is the momentum operator with com- 
mutator given by [x, p] = i. Moreover, it is obvious 
from the Hamiltonian that the following commutation 
relations hold [H, x\ ~ —ip, [H, p\ — ix. Identifying 
Q and P as x and p, we see that 



[H, (x, p)] = (x, p) 



i 
-i 



(5) 



Diagonalizing the matrix on the right hand side of 
eq(|) gives 

[H, {x, p)U] = (x, p)UU-' (^^q}U 



so that 



1 /I 1 



Thus, 



[H, (a, at)] = (a, a^) ( J ) , 



(6) 



III. SHIFT OPERATORS 

We next see how the shift operators can be ob- 
tained for Hamiltonians with arbitrary solvable poten- 
tials. This method can also be regarded as a means 
of searching systematically for the solvable potentials. 
To this end, we begin with the following Hamiltonian 



(7) 



where X(x) is an arbitrary function of position and 
V{x) is an arbitrary potential. Define 



P{x,p) = Yix)— + Zix) 
ax 



(8) 



with arbitrary functions Y{x) and Z{x) to be deter- 
mined. It is easy to show that the commutation rela- 
tion between H and P is 

[H, P] = {2X{x)Y\x) - X\x)Y{x))^ 

+X{x){Y"{x) + 2Z'{x))4- + X{x)Z"{x) 
ax 



Y{x)V'{x). 



(9) 



and we can therefore identify a — "^/^(^ ^ ^P) ^'^'^ 

= --i=(x — ip) as the natural shift operators aris- 
ing from the diagonalization process. The harmonic Furthermore, if we impose for simplicity 



Since X{x), Y{x) and Z{x) are arbitrary functions, it 
is possible to set 

X{x){Y"{x) + 2Z'{x)) = aY{x) (10) 
2X{x)Y'{x) - X'{x)Y{x) = {(3Q{x) + j)X{x), (11) 

where a, (3 and 7 are arbitrary constants and Q{x) 
is the generalized "coordinate" to be determined, so 
that 

[H, P] = {(iQ{x) + j)H + aP- {f3Q{x) + -f)V{x) 

-aZ{x) ~ Y{x)V'{x) + X{x)Z"{x). (12) 
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X{x)Z"{x) - if3Q{x)+j)V{x) - aZ{x) 

-Y{x)V'{x)^Q{x), (13) 

and solved for the operator Q{x) as 



- pVix) 



{X{x)Z"{x) - -yVix) 

- aZ{x) -'Y{x)V'ix)), (14) 



we find 



[H, P] = Q{x){f3H + l) + aP + -fH. (15) 

Concomitantly, we can also work out explicitly the 
expression for the commutator of H with Q{x) giving 



[H, Q{x)] = 2X{x)Q'ix)— + X{x)Q"{x). (16) 
ax 



Finally, one can set 

X{x)Q'{x) = \Y{x), 

so that 

[H, Q{x)] = 2XP - 2XZ{x) + X{x)Q"ix). 
If the further constraint 

- 2XZ{x) + X{x)Q"{x) = iyQ{x) + r, 
is imposed, we arrive at 

[H, Q{x)] = 2\P + vQ{x)+T. 



(17) 



(18) 



(19) 



Here, A, v and r are constants to be determined. Note 
that the quantities r and can be absorbed in cer- 
tain cases by a redefinition of the operators Q{x) and 
P so that we finally obtain a semi-closed algebra. In 
the process, the parameters a, 7, A, v, t are not 
arbitrary but must satisfy some constraints. It is also 
interesting to note that in the above construction of 
the shift operators, we are not constrained by com- 
mutation relation between Q(x) and P . Moreover, we 
do not need to require the operators iJ, Q and P form 
a closed algebra, although for most of the well known 
exactly solvable Hamiltonians, these operators always 
inevitably form a closed algebra. However, to obtain 
the ground state, we need to require that the Hamil- 
tonian factorizes according to 5*1 6*2 ~ F{H, li, Ci) for 
some arbitrary function, F. 

In the construction, we should bear in mind that 
Y{x), Z{x),Q{x) and a, 7, A, t are all de- 
termined by X(x) and V{x) when we construct the 
shift operators for a given Hamiltonian. Conversely, 
if we proceed with a systematic search for exactly solv- 
able potentials based on this method, we can have as 
many degrees of freedom as there are free parameters 
and functions, namely X{x), Y(x), Z{x),Q{x) and 
a, 1, A, I/, T. 



IV. SOME SOLVABLE POTENTIALS 

We can perform a systematic search for a class of 
exactly solvable potentials by solving the equations 
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,|1§|). Furthermore, the method is tanta- 
mount to a construction the corresponding shift oper- 
ators so that it is possible to obtain the energy spec- 
trum and the underlying algebraic formulation if nec- 
essary. Indeed, most of well known solvable potentials 
can be reproduced using the method by an appropri- 
ate choice of functions and parameters. 
Case 1. X{x) = -1, Y{x) = 1. 
On substituting into the equations (|l^Jll|,|l|,|l|,^ , 
consistency requires the imposition of the following 
constraints 

Q'{x) = -A, 
-PQix) -7 = 0, 
-2Z'(x) = a, 
-2XZ{x) - Q"{x) = uQ{x) + r, 
V'{x) + {pQix) + 7)F(x) = ~Z"{x) - aZ{x) ~ Q{x). 

It is easy to solve these constrained equations to ob- 
tain 



Q{x) = -Xx + ci, Z{x) = -—x + C2, 



and (3 = 0, 7 = 0, v = —a and t = civ — 2Xci 
with the integral constants Ci,i = 1,2,3. This is the 
harmonic oscillator. 

By substituting the solution of the constrained 
equations into eq(^), we get 

[i?, Q] = -olQ + 2AP - (2Ac2 -f vex), 
[H, P] = Q + aP. 

A simple redefinition of the operators Q and P using 

Q = Q + f and P = P + g where / = — a^^^— ^--77— 

a'^ + 2X 

and g = — '^^ ~^ so that the operators Q and P 
-I- 2A 

take the form 

Q = -Xx + ci+ a — „ 

a-^ -I- zA 

~ d a 2Xc2 + vci 

~ 2 ~ a2 + 2A 

immediately give the equations 

[H, Q] = -aQ + 2XP, 
[H, P]=Q + aP, 
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or more succinctly as 



-a 1 
2A a 



which can be easily diagonalized to give the shift op- 
erators Si and 5*2 which satisfies the commutation re- 
lation 



[H,{Si,S2)]^{Si,S2) 
with 

Si - ■ 



2A 







\/a2 + 2A 



1 



Va2 + 2A 
1 



-f + 2A 



Q 



p 
p 



and 



Va2 + 2A 



(20) 



It is clear from eq(pO|) that the shift operators 5*1 and 
S2 are really the usual lowering and raising operators 
resepctively for the harmonic oscillators ( simply set 



A = 



and = 



-). Hav- 



(a2 + 2A)'/^ (a2 + 2A)'/^' 

ing constructed the shift operators, the ground state 
eigenvalue and eigenvector can be solved by the ac- 
tion of the lowering operator on a lowest weight state 
using the equations 



The results are easily found to be 

^i?ix2 



tpo = B3 exp ^ 



Eo 



2A + C3 - 



B2X 
(ci + C2a)2 



2A 



(21) 

(22) 
(23) 



where Bi, i = 1, 2, 3 are constants of integration. 
Case 2. X{x) = -1, Y{x) = x. 
As in Case 1, a simple substitution into the equa- 



tions (|l^,|ll|,|l^,^ 18) requires for consistency the fol- 
lowing constraints 

Q' {x) = —Ax, 

-/3Q(a:) - 7 = -2, 
—2Z'{x) = ax, 
-2XZ{x) - Q"{x) = iyQ{x) + r, 
xV'{x) + iPQ{x) + -f)V{x) = -Z"{x) - aZ{x) - Q{x) 

By solving the equations, we obtain 



Qix) 



A 



x'^ + ci, Z{x) 



C2, 



with /3 = 0, 7 = 2, V = —a and r = A(l — 2c2) + aci. 
This gives us the radial harmonic oscillator potential. 
The rest of the details are similar to Case 1, but we 
will provide some of the crucial results. 

The shift operators Si and 5*2 can be found through 
a similar method with the explicit form of the opera- 
tors given by 



51 = Q + (a - Vq^^ + 2A)P 

52 = g + (a + v^a^ + 2A)P 

as before but with the redefined operators Q and P 
defined by 



Q 



A 2 

'—X 

2 



4A 



ci 



2A 

a 



-H 



P = X- x^ 

dx 4 



:e2 -I-2A 
2a 



C2 



(A + aci - 2AC2) 
H 



2-I-2A" 

(X + aci -2AC2). 



2A' 



Moreover, the following commutation relation, similar 
to Case 1, is also satisfied 



[H,{Si,S2)] = {Si,S2) 



-^/a^ + 2X 








Va2 + 2A 



Note that these operators depend on H so that the 
eigenvalues associated with the raising and lowering 
operators depend on the energy level it is acting on. 
Thus the tower of states is no longer equally spaced. 
The commutation relation for the shift operators is 

[Si,S2] = -^4L=i2H + CI+ C2a - ^). (24) 
Va + 2A ^ 

The ground state can be shown to be 

V'o = i?4a;-(^^^°+^^) eM~^x^) (25) 
with energy given implicitly by 

{B2Eo + B3f + {B2Eo + B3)-C3 = 0. (26) 



Case 3. X{x) = -1, Y{x) = ae""" + foe"^^. 
A similar substitution into equations 



( |lO| , |ll|jl4yi7| , |lS| ) and solving the resulting consistency 
equations give the solution 



Q{x) 



A 



(ae'=^-6e-'=^)-|-ci. 



n-) = 3(A+2)V + ^ + -(--«C2-cO 



Z{x) 



2c 



-iae'^ + be-'^) + C2, 
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V{x) = C3(ae=" + 6e-=")-2 + —^[{a + c^f + 2A] 

and /? = -fc^, 7 = ^ci, = -a - 2c^, t = 
— 2Ac2 — i^ci. This is the generahzed Poschl- Teller po- 
tential Moreover, we also have the the relation 
ci + ac2 = 0. 

The shift operators can be found similarly using the 
technique described in Case 1 and the commutation 
relations for the shift operators, 5*1 and 5*2, are 

{H, 5i] = 5i (-c^ - V(a + c2)2 + 2A-4c2ij) (27) 

[ff, 52] - 52 (-c2 + V(a + c2)2 + 2A - 4c2i/) . (28) 

The explicit form of the shift operators are given by 

51 = -Q7^((a + c2) - V(a + c2)2_4c2i/ + 2A) + P 

52 = -Q7^((a + c') + V(« + c2)2-4c27J + 2A) + P 



where 



P = P - C2. 

with commutation relation 



[Q, P] = A (ae'^^ + foe" 



and 



[5i,52] = -8a6cy (^±^)2 + ^ - (29) 
The ground state is given by 



V>o = C4(ae^^ + &e-^-) V '''' (30) 

and the ground state energy Pq can be determined 
through the equation 



a + c2 A 



, , a + c2 ^ „ A C4 , ^ 

-V(^)^+2;^-^°+4i=°- 

Case 4. X(a;) = —1, Y{x) = asin(fca;) + 6cos(fcx). 
For this case, we get 

Q(x) — —(acos(kx) — 6sin(fcx)) + Ci, 
k 



and /? = 7 = -|cifc2^ = 2fc2 _ r = 

ci{a — 2fc2) — 2Ac2. This is the familiar Poschl- Teller 
potential. 

The shift operators are given by 



[P, Si] = 5i(fc2 - v/(a-fc2)2 + 2A-F4fc2i?) 
[P, 52] = 52(fc2 + v/(a-fc2)2 + 2A + 4fc2p) 

with the operators 5i and 52 defined by 
5i = (asin(fca;) -f 6cos(fca;)) 



dx 



1 

2k 



(acos(fcx) - &sin(fca;))^(a - fc2)2 + 2A + 4fc2p 



52 — (asin(fca;) + bcos(kx)) 



1 



dx 



+ Trr{acos(kx) - fosin(fca::)) ^(a - fc2)2 + 2A + 4fc2p 
2k 

and the commutator given by 
[5i,52] = -2fc(a2 + 62 



The ground state, -00, can be solved and found to be 
given by the formula 



4>o = C4(asin(fca;) -|-6cos(fca;))ii^V("-fc')'+2^+4fe^£^o 

(33) 

where Pq is the ground state energy determined 
through the equation 



4fc2 



62 



C3 + 2A:V(" - k^)^ + 2A + 4fc2Po 

- ((a- fc2)2 + 2A-|-4fc2£;o) =0. (34) 



Case 5. X{x) = —x, Y{x) = x. 

A similar substitution of above trial functions yields 



Q{x) = — Aa; + ci, Z{x) 



:x + C2, 



y{x) = \{X + ^)x + — ~{ci+ ac2) 
L L X 

and /? = 0, 7 = 1, v ^ —a, r — ac\ — 2Ac2. This is 
the case of Coulomb potential for the hydrogen atom. 

The shift operators, S\ and 52, satisfy the commu- 
tation relations 



[P, 5i] = -Va^ + 2A5i 
[P,52] = Va2-f2A52, 



Zix) = — — — (acos(fca;) — 6sin(fcx)) +C2, 
2k 



V{x) — C3{asm{kx) + bcos{kx)) 2 — 



(fc2 - a)2 + 2A 
4fc2 



in the same manner as the case for the harmonic os- 
cillator in Case 1 and the radial harmonic oscillator 
in Case 2 but with the operators defined as 

(5i,52) = (Q,P)[/ 
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where 



A X 2A Ci + ac2 

Q = -^^+ 2 , ox -^ + 2A ^ ^ 



5 d a 
F = X— —X 



a ci + ac2 

-H + a- 



dx 2 a2 + 2A a2 + 2A ' 

and the matrix U given by 

U = ( a-Va2+2A a+Va2+2A ^ 

It is instructive to note that the commutation relation 
for the shift operators can be computed succinctly as 



+ 2A 



{H + ci+ ac2). 



A direct computation using SiipQ = yields the 
ground state as 



where the ground state energy Eq is given by 



(35) 



£;o = ^^^^^^(l±\/m^)-ci-ac2, (36) 

choosing appropriately the correct minimal value. 

Case 6. X{x) = -e^^, Y{x) = 1. 

As a final case, wc look the above choice, so that 
upon substitution into the appropriate equations gives 

Q{x) = -e-""" + ci, Z{x) = ^e-""" + C2, 
c 2c 



V{x) 



2A ' -2 



_i!:e-c^ + e3e- + |(2c2-c) + £l 



and/3 = 0, 7 = — c, = —a, r = aci — A(c+2c2). We 
note that this is simply the case of Morse potential. 
As in Case 5, the shift operators, and 52, obey the 
relations 

[H,Si] = -\/a2 + 2A5'i 
[H, S2] = ^o? + 2A52, 

with the shift operators defined by 

{8^,82) = {Q,P)U 

where 

c + 2A 

da;^2c a2 + 2A ' 

and the matrix U given by 

1 1 

TJ = ( a-Va^+2A a+Va2+2A 
1 1 



as before. The commutation relation between the shift 
operators, >Si and 6*2, now reads 



[SuS2 



2c2 ^(^^2c,+a(c + 2c2)^_ ^3^^ 



\/a2 _^ 2A 



2c 



Moreover, a straightforward calculation yields the 
ground state as 



V'o = C4 exp 



{- 



Vq!^ _^ 2A 
2c2 ' 



, a{c + 2c2] + cVa2 + 2A + 2ci - 2ci;o , 
+ ^ _ . „ X ) , (38) 



2Va2 + 2A 

and the ground state energy, Eq, is determined by 

Eq = ^{2ci +a(c + 2c2) 

+c\/a2 + 2A ± V2c3(a2+2A)} , (39) 

with an appropriate choice of sign. 

Finally, we would like to comment briefly on the 
last two cases. Normally, for radial Coulomb prob- 
lem and Morse oscillator, we deal with the following 
Schrodinger equations 



+ ^)V'^,i(r-) =0 



r 



and 



(-^ + iHl - e-yf - = -{l-v- i)2^,,, 

respectively. However, it is always possible to trans- 
form the above equations into more convenient forms 
as 



with p = —r and 



^2 l{l + 1) 



(-e"^ + {l-v- i)2e- + e--)V.,; = 2ZV.,;. 

To get a deciper insight into this transformation, we 
note that the Hamiltonian concerned can be rewritten 
as 

where T{x) and J?(a;) are suitably chosen functions. 
The eigenvalue problem for stationary eigenstates 
takes the form 

(-^ + R{x) + = E^(x). (40) 

Multiplying by the function T(x), we see that the 
same eigenvalue problem now appear as 
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+ Tix){R{x) - E)Wx) = -^P{x). (41) 

Hence, the original eigenvalue problem is now trans- 
formed into an equivalent eigenvalue problem with the 
Hamiltonian 



H' = -T{x)— + T{x){R{x) - E). 



dx 



(42) 



It is precisely with the new form that the potential 
can be solved more easily. For Coulomb and Morse 
potentials, we have therefore used the modified Hamil- 
tonian. 



V. DISCUSSION AND CONCLUSION 

It is known that the local behavior of most solvable 
potentials reduces either to the harmonic oscillator or 
the Poschl- Teller potential |^^. This fact is further 
confirmed by the commutation relations between the 

^ A 



Hamiltonian and the shift operators in section IV 



quick glance at the various cases clearly shows that 
the general form for the commutation relations can 
be classified into two major categories, those belong- 
ing to the harmonic oscillator (Case 1, Case 2, Case 
5 and Case 6), corresponding to /3 = 0, and those 
belonging to the Poschl- Teller potential (Case 3 and 
Case 4), corresponding to /? ^ 0. Thus, algebraically, 
we should expect similar behavior for potentials be- 
longing to the same class. 

The inherent nonlinearity in the algebraic struc- 
ture of the shift operators, in which Fi{H)Si — 
SiGi{H),i ~ 1,2, should have a deep physical con- 
nections to the study of exactly solvable potentials in 
supersymmetric theory [|l3|, Moreover, the tech- 
nique presented can be systematically analyzed to en- 
compass some of these supersymmetric potentials. In- 
deed, in section IV, we have not entirely exhausted all 
possibilities and provided a comprehensive list. How- 
ever, it is in principle possible to do so. By exploring 
more complicated functions X{x), Y(x) and Z(x), it 
is expected that more complicated solvable potentials 
may arise. 

Recently, there has been much interest generated 
regarding the computation of Franck-Condon factors 
for anharmonic oscillators j2^, In particular the 
Morse and Poschl- Teller potentials do not in general 
provide a simple analytical expression for the Franck- 
Condon factors [pi , 32 1 . Our technique should be able 



to cast new light on the exact computation of the 
Franck-Condon factors for such potentials. More work 
in this direction will be reported elsewhere. Finally, 
we note that the method can be extended easily to 
many-body potentials. This work was supported pri- 
marily by NUS Research Grant No. RP3982713. In 
addition, we would also like to thank the referee for 
his invaluable comments. 
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